Impact of genotype-calling methodologies on genome-wide association and genomic prediction in polyploids.

Discovery and analysis of genetic variants underlying agriculturally important traits are key to molecular breeding of crops. Reduced representation approaches have provided cost-efficient genotyping using next-generation sequencing. However, accurate genotype calling from next-generation sequencing data is challenging, particularly in polyploid species due to their genome complexity. Recently developed Bayesian statistical methods implemented in available software packages, polyRAD, EBG, and updog, incorporate error rates and population parameters to accurately estimate allelic dosage across any ploidy. We used empirical and simulated data to evaluate the three Bayesian algorithms and demonstrated their impact on the power of genome-wide association study (GWAS) analysis and the accuracy of genomic prediction. We further incorporated uncertainty in allelic dosage estimation by testing continuous genotype calls and comparing their performance to discrete genotypes in GWAS and genomic prediction. We tested the genotype-calling methods using data from two autotetraploid species, Miscanthus sacchariflorus and Vaccinium corymbosum, and performed GWAS and genomic prediction. In the empirical study, the tested Bayesian genotype-calling algorithms differed in their downstream effects on GWAS and genomic prediction, with some showing advantages over others. Through subsequent simulation studies, we observed that at low read depth, polyRAD was advantageous in its effect on GWAS power and limit of false positives. Additionally, we found that continuous genotypes increased the accuracy of genomic prediction, by reducing genotyping error, particularly at low sequencing depth. Our results indicate that by using the Bayesian algorithm implemented in polyRAD and continuous genotypes, we can accurately and cost-efficiently implement GWAS and genomic prediction in polyploid crops.

Effects of tall alleles SD1-in and SD1-ja to the dwarfing allele sd1-d originating from 'Dee-geo-woo-gen' on yield and related traits on the genetic background of indica IR36 in rice.

sd1-d originating from 'Dee-geo-woo-gen' has been utilized to develop short-culmed indica varieties adaptable to higher fertilizer-application. Its tall alleles SD1-in and SD1-ja are harbored in indica and japonica subspecies, respectively. The sd1-d of indica IR36 was substituted with SD1-in or SD1-ja by recurrent backcrossing with IR36, and two tall isogenic lines ("5867-36" and "Koshi-36") were developed. IR36, 5867-36 and Koshi-36 were grown in a paddy field in three years, and yield and related traits were measured, the effects of SD1-in and SD1-ja on yielding ability and related characteristics were examined on the genetic background of IR 36. SD1-in decreased panicle number per m2 but increased spikelet number per panicle, ripened-grain percentage and 1000-grain weight, compared with sd1-d, resulting in the increase of yield. The increase of 1000-grain weight by SD1-in, caused by the increases of length, width and thickness of grain, was due to the increases of the length and width of lemma. SD1-ja did not significantly affect yield, mainly because the decrease of panicle number per m2 was compensated by the enlarged 1000-grain weight owing to the increase of lemma length. Serious lodging was observed in long-culmed 5867-36, suggesting that sd1-d is indispensable for indica breeding programs.

Characterization of the Ghd8 Flowering Time Gene in a Mini-Core Collection of Miscanthus sinensis.

The optimal flowering time for bioenergy crop Miscanthus is essential for environmental adaptability and biomass accumulation. However, little is known about how genes controlling flowering in other grasses contribute to flowering regulation in Miscanthus. Here, we report on the sequence characterization and gene expression of Miscanthus sinensisGhd8, a transcription factor encoding a HAP3/NF-YB DNA-binding domain, which has been identified as a major quantitative trait locus in rice, with pleiotropic effects on grain yield, heading date and plant height. In M. sinensis, we identified two homoeologous loci, MsiGhd8A located on chromosome 13 and MsiGhd8B on chromosome 7, with one on each of this paleo-allotetraploid species' subgenomes. A total of 46 alleles and 28 predicted protein sequence types were identified in 12 wild-collected accessions. Several variants of MsiGhd8 showed a geographic and latitudinal distribution. Quantitative real-time PCR revealed that MsiGhd8 expressed under both long days and short days, and MsiGhd8B showed a significantly higher expression than MsiGhd8A. The comparison between flowering time and gene expression indicated that MsiGhd8B affected flowering time in response to day length for some accessions. This study provides insight into the conserved function of Ghd8 in the Poaceae, and is an important initial step in elucidating the flowering regulatory network of Miscanthus.

Managing flowering time in Miscanthus and sugarcane to facilitate intra- and intergeneric crosses.

Miscanthus is a close relative of Saccharum and a potentially valuable genetic resource for improving sugarcane. Differences in flowering time within and between Miscanthus and Saccharum hinders intra- and interspecific hybridizations. A series of greenhouse experiments were conducted over three years to determine how to synchronize flowering time of Saccharum and Miscanthus genotypes. We found that day length was an important factor influencing when Miscanthus and Saccharum flowered. Sugarcane could be induced to flower in a central Illinois greenhouse using supplemental lighting to reduce the rate at which days shortened during the autumn and winter to 1 min d-1, which allowed us to synchronize the flowering of some sugarcane genotypes with Miscanthus genotypes primarily from low latitudes. In a complementary growth chamber experiment, we evaluated 33 Miscanthus genotypes, including 28 M. sinensis, 2 M. floridulus, and 3 M. ×giganteus collected from 20.9° S to 44.9° N for response to three day lengths (10 h, 12.5 h, and 15 h). High latitude-adapted M. sinensis flowered mainly under 15 h days, but unexpectedly, short days resulted in short, stocky plants that did not flower; in some cases, flag leaves developed under short days but heading did not occur. In contrast, for M. sinensis and M. floridulus from low latitudes, shorter day lengths typically resulted in earlier flowering, and for some low latitude genotypes, 15 h days resulted in no flowering. However, the highest ratio of reproductive shoots to total number of culms was typically observed for 12.5 h or 15 h days. Latitude of origin was significantly associated with culm length, and the shorter the days, the stronger the relationship. Nearly all entries achieved maximal culm length under the 15 h treatment, but the nearer to the equator an accession originated, the less of a difference in culm length between the short-day treatments and the 15 h day treatment. Under short days, short culms for high-latitude accessions was achieved by different physiological mechanisms for M. sinensis genetic groups from the mainland in comparison to those from Japan; for mainland accessions, the mechanism was reduced internode length, whereas for Japanese accessions the phyllochron under short days was greater than under long days. Thus, for M. sinensis, short days typically hastened floral induction, consistent with the expectations for a facultative short-day plant. However, for high latitude accessions of M. sinensis, days less than 12.5 h also signaled that plants should prepare for winter by producing many short culms with limited elongation and development; moreover, this response was also epistatic to flowering. Thus, to flower M. sinensis that originates from high latitudes synchronously with sugarcane, the former needs day lengths >12.5 h (perhaps as high as 15 h), whereas that the latter needs day lengths <12.5 h.

Training Population Optimization for Genomic Selection in Miscanthus.

Miscanthus is a perennial grass with potential for lignocellulosic ethanol production. To ensure its utility for this purpose, breeding efforts should focus on increasing genetic diversity of the nothospecies Miscanthus × giganteus (M×g) beyond the single clone used in many programs. Germplasm from the corresponding parental species M. sinensis (Msi) and M. sacchariflorus (Msa) could theoretically be used as training sets for genomic prediction of M×g clones with optimal genomic estimated breeding values for biofuel traits. To this end, we first showed that subpopulation structure makes a substantial contribution to the genomic selection (GS) prediction accuracies within a 538-member diversity panel of predominately Msi individuals and a 598-member diversity panels of Msa individuals. We then assessed the ability of these two diversity panels to train GS models that predict breeding values in an interspecific diploid 216-member M×g F2 panel. Low and negative prediction accuracies were observed when various subsets of the two diversity panels were used to train these GS models. To overcome the drawback of having only one interspecific M×g F2 panel available, we also evaluated prediction accuracies for traits simulated in 50 simulated interspecific M×g F2 panels derived from different sets of Msi and diploid Msa parents. The results revealed that genetic architectures with common causal mutations across Msi and Msa yielded the highest prediction accuracies. Ultimately, these results suggest that the ideal training set should contain the same causal mutations segregating within interspecific M×g populations, and thus efforts should be undertaken to ensure that individuals in the training and validation sets are as closely related as possible.

Population structure of Miscanthus sacchariflorus reveals two major polyploidization events, tetraploid-mediated unidirectional introgression from diploid M. sinensis, and diversity centred around the Yellow Sea.

BACKGROUND AND AIMS: Miscanthus, a C4 perennial grass native to East Asia, is a promising biomass crop. Miscanthus sacchariflorus has a broad geographic range, is used to produce paper in China and is one of the parents (along with Miscanthus sinensis) of the important biomass species Miscanthus × giganteus. The largest study of M. sacchariflorus population genetics to date is reported here. METHODS: Collections included 764 individuals across East Asia. Samples were genotyped with 34 605 single nucleotide polymorphisms (SNPs) derived from restriction site-associated DNA sequencing (RAD-seq) and ten plastid microsatellites, and were subjected to ploidy analysis by flow cytometry. KEY RESULTS: Six major genetic groups within M. sacchariflorus were identified using SNP data: three diploid groups, comprising Yangtze (M. sacchariflorus ssp. lutarioriparius), N China and Korea/NE China/Russia; and three tetraploid groups, comprising N China/Korea/Russia, S Japan and N Japan. Miscanthus sacchariflorus ssp. lutarioriparius was derived from the N China group, with a substantial bottleneck. Japanese and mainland tetraploids originated from independent polyploidization events. Hybrids between diploid M. sacchariflorus and M. sinensis were identified in Korea, but without introgression into either parent species. In contrast, tetraploid M. sacchariflorus in southern Japan and Korea exhibited substantial hybridization and introgression with local diploid M. sinensis. CONCLUSIONS: Genetic data indicated that the land now under the Yellow Sea was a centre of diversity for M. sacchariflorus during the last glacial maximum, followed by a series of migrations as the climate became warmer and wetter. Overall, M. sacchariflorus has greater genetic diversity than M. sinensis, suggesting that breeding and selection within M. sacchariflorus will be important for the development of improved M. × giganteus. Ornamental M. sacchariflorus genotypes in Europe and North America represent a very narrow portion of the species' genetic diversity, and thus do not well represent the species as a whole.

Contrasting allelic distribution of CO/Hd1 homologues in Miscanthus sinensis from the East Asian mainland and the Japanese archipelago.

The genus Miscanthus is a perennial C4 grass native to eastern Asia and is a promising candidate bioenergy crop for cool temperate areas. Flowering time is a crucial factor governing regional and seasonal adaptation; in addition, it is also a key target trait for extending the vegetative phase to improve biomass potential. Homologues of CONSTANS (CO)/Heading date 1(Hd1) were cloned from Miscanthus sinensis and named MsiHd1. Sequences of MsiHd1 homologues were compared among 24 wild M. sinensis accessions from Japan, 14 from China, and three from South Korea. Two to five MsiHd1 alleles in each accession were identified, suggesting that MsiHd1 consists of at least three loci in the Miscanthus genome. Verifying the open reading frame in MsiHd1, they were classified as putative functional alleles without mutations or non-functional alleles caused by indels. The Neighbor-Joining tree indicated that one of the multiple MsiHd1 loci is a pseudogene locus without any functional alleles. The pseudogene locus was named MsiHd1b, and the other loci were considered to be part of the MsiHd1a multi-locus family. Interestingly, in most Japanese accessions 50% or more of the MsiHd1a alleles were non-functional, whereas accessions from the East Asian mainland harboured only functional alleles. Five novel miniature inverted transposable elements (MITEs) (MsiMITE1-MsiMITE5) were observed in MsiHd1a/b. MsiMITE1, detected in exon 1 of MsiHd1a, was only observed in Japanese accessions and its revertant alleles derived from retransposition were predominantly in Chinese accessions. These differences in MsiHd1a show that the dependency on functional MsiHd1a alleles is different between accessions from the East Asian mainland and Japan.

Is RNA editing implicated in group II intron survival in the angiosperm mitochondrial genome?

Introns may be considered as optional because they are removed from mRNA molecules, but introns are fairly preserved for unknown reasons. Previously, the mitochondrial rps3 gene of sugar beet ( Beta vulgaris L., Caryophyllales) was shown to represent a unique example of an intron loss. We have determined the distribution of the rps3 intron in 19 Caryophyllalean species. The intron was absent from the Amaranthaceae and the Achatocarpaceae. In the Caryophyllaceae, Dianthus japonicus rps3 was pseudogenized, but the intronic sequence was retained. Intact intron-bearing rps3 copies were cloned from Portulaca grandiflora and Myrtillocactus geometrizans , members of the sister clade of the Amaranthaceae-Achatocarpaceae-Caryophyllaceae clade. Most of the C-to-U RNA-editing sites in P. grandiflora and M. geometrizans rps3 transcripts were homologous in the two species, as well as in the sugar beet rps3, which, unlike the other 12 rps3 transcripts, lacks editing in the exonic regions around the intron. Provided that the loss of editing preceded the loss of rps3 intron, it appears conceivable that a requirement for editing could have prevented the loss of group II introns retained in angiosperm mitochondrial genomes. This interpretation is an alternative to the conventional one that views the loss of editing as a mere trace of RNA-mediated gene conversion.

Differentiation in wild-type allele of the sd1 locus concerning culm length between indica and japonica subspecies of Oryza sativa (rice).

A dwarfing gene (allele) sd1-d has been intensively utilized to develop short-culm indica varieties in southeast Asia up to now. Before the first sd1-d-carrying variety IR8 was released, rice researchers had recognized the general tendency that culm length is higher in indica varieties than in temperate-japonica ones. Inter-subspecific difference of the tall (wild-type) allele SD1 at the sd1 locus was examined on the common genetic background, using five isogenic lines developed by substituting sd1-d of the recurrent parent IR36 by SD1s of two indica varieties, two temperate-japonica varieties and one tropical-japonica variety. The two indica -donor isogenic lines had longer culms than the three japonica-donor isogenic lines consistently in two different environmental conditions. Moreover, nonsynonymous single-nucleotide polymorphism between the two subspecies was detected at two sites in Exon 1 and Exon 3 of the sd1 locus. It is demonstrated that the inter-subspecific differentiation of SD1 contributes height difference between indica and japonica. The indica-originating and japonica-originating alleles at the sd1 locus were designated as SD1-in(t) and SD1-ja(t), respectively.

Large 3' UTR of sugar beet rps3 is truncated in cytoplasmic male-sterile mitochondria.

Genomic alteration near or within mitochondrial gene is often associated with cytoplasmic male sterility (CMS). Its influence on the expression of the mitochondrial gene was proposed as one of the possible causes of CMS. In sugar beet mitochondrial rps3, whose downstream 1,056-bp region contains Norf246, an apparently non-functional open reading frame (ORF), was deleted in CMS mitochondria. In our previous study, normal rps3 (3.8 kb), CMS rps3 (2.7 kb), and Norf246 (3.8 and 0.9 kb) were shown to be transcribed. The present study was conducted to determine whether the deletion affected gene expression. Reverse transcription (RT)-PCR analysis revealed the co-transcription of rps3 and Norf246. By circularized RNA (CR) RT-PCR analysis, the 5' and 3' termini of the 3.8- and the 0.9-kb transcripts were determined. The results suggested that the 3.8-kb transcripts were the rps3 mRNA bearing ~464-base 5' untranslated region (UTR) and ~1,508-base 3' UTR, whereas no functional ORF was observed in the 0.9-kb transcripts. CR-RT-PCR revealed that the 3' UTR of the 2.7-kb transcripts was reduced to ~460 bases. However, no difference in the accumulation of RPS3 polypeptide and RNA editing was detected by protein gel blot analysis and cDNA sequencing. Although the deleted region encoded the truncated-atp9 that was edited, no influence on the pattern and frequency of RNA editing of genuine atp9 was evident. The results eliminated rps3 as a candidate for the CMS gene, making preSatp6, a unique ORF fused with CMS atp6, the sole CMS-associated region in sugar beet.

Diversification in flowering time due to tandem FT-like gene duplication, generating novel Mendelian factors in wild and cultivated rice.

The complex structure of a single Mendelian factor widespread in the Asian cultivated rice (Oryza sativa) and its wild progenitor (Oryza rufipogon) that caused diverse phenotypes in the timing of flowering under natural field conditions was investigated in near isogenic lines. These near isogenic lines showed differences in flowering time despite all eight accessions collected from tropical regions possessing a recessive gene allelic to the se-pat gene. Fine mapping in two of these near-isogenic lines revealed that cultivated (Patpaku) and wild (W593) accessions had three and two linked quantitative trait loci (QTL) in the candidate regions, respectively, showing that Patpaku and W593 possessed linked QTLs with different effects in addition to the commonly-observed recessive gene (se-pat). Molecular dissection suggested that the tandemly duplicated FT-like genes (Hd3a and RFT1) could be the candidate genes for these QTLs. Interestingly, the linked QTLs differed in their epistases, degree of dominance, and genotype x environment interactions. The nucleotide sequences showed that RFT1 has diverged more rapidly than Hd3a during rice evolution, suggesting phenotypic diversification of the two genes. Phylogenetic analysis implied that the se-pat(+) alleles might have emerged in different lineages within O. sativa. The present results strongly suggest that nucleotide divergence and shuffling of the linked QTLs by recombination might have created novel Mendelian factors that probably contribute to responding to local environments.

Evidence for an evolutionary force that prevents epigenetic silencing between tail-to-tail rice genes with a short spacer.

During the course of evolution, the genome should have toned down various types of genomic noise, such as those that cause the unstable expression or gene silencing observed in transgenic organisms. We found a rice genomic segment where two genes, encoding 5-enolpyruvylshikimate-3-phosphate synthase (EPSPs) and ribosomal protein small subunit 20 (rps20), are located in a tail-to-tail orientation and separated by only 300 bp of spacer. It is possible that this kind of structure would give rise to unstable expression due to antisense RNA derived from the neighboring gene. We examined this possibility using Northern blot, reverse transcription-polymerase chain reaction (RT-PCR), and 3' RACE analyses, but obtained no evidence for instability or antisense RNAs of these housekeeping genes. Comparison of the sequences in the corresponding regions among related rice species revealed a lower level of genetic divergence of both the 3'-untranslated region (3'-UTRs) than of the other noncoding regions; in particular both of the boundaries between the 3'-UTRs and the spacer were markedly conserved. The conservation of both the terminal regions is most likely the result of purifying selection, implying a functional role for the strict termination of the transcription of these genes to prevent gene-silencing-related events.

Genealogy of the "Green Revolution" gene in rice.

During the "Green Revolution" of rice, high-yielding varieties (HYVs) were developed using a semi-dwarf gene (sd1 or OsGA20ox2). The presence or absence of the two mutant alleles (DGWG type in Dee-geo-woo-gen and JKK type in Jikkoku) were surveyed by PCR using 256 accessions of eight wild and two cultivate rice species. The DGWG allele was detected in a landrace (Oryza sativa) and two accessions of wild rice (O. rufipogon), all of which are from China, showing their limited distribution. Genealogical studies of the OsGA20ox2 gene showed that the 62 sequences of O. sativa and O. rufipogon included 20 distinct haplotypes, indicating that the species complex contained OsGA20ox2 genes from two different lineages. The silent site nucleotide diversities (pi and theta(w)) were extremely low in Japonica rice, suggesting a genetic bottleneck. The haplotype network showed that the DGWG and JKK alleles were derived in different lineages. The DGWG carrier (W1944) had unique polymorphisms in the surrounding region of the locus, suggesting that the DGWG allele has been preserved in the wild progenitor, rather than that the DGWG allele has been introgressed from HYVs to W1944. Although a semi-dwarfing plant is a weak competitor under saturated fields, the crossing experiment revealed that the DGWG variant might have been preserved as a hidden variation in the genetic background of wild rice, without expressing a short-stature.

Reconstruction of putative DNA virus from endogenous rice tungro bacilliform virus-like sequences in the rice genome: implications for integration and evolution.

BACKGROUND: Plant genomes contain various kinds of repetitive sequences such as transposable elements, microsatellites, tandem repeats and virus-like sequences. Most of them, with the exception of virus-like sequences, do not allow us to trace their origins nor to follow the process of their integration into the host genome. Recent discoveries of virus-like sequences in plant genomes led us to set the objective of elucidating the origin of the repetitive sequences. Endogenous rice tungro bacilliform virus (RTBV)-like sequences (ERTBVs) have been found throughout the rice genome. Here, we reconstructed putative virus structures from RTBV-like sequences in the rice genome and characterized to understand evolutionary implication, integration manner and involvements of endogenous virus segments in the corresponding disease response. RESULTS: We have collected ERTBVs from the rice genomes. They contain rearranged structures and no intact ORFs. The identified ERTBV segments were shown to be phylogenetically divided into three clusters. For each phylogenetic cluster, we were able to make a consensus alignment for a circular virus-like structure carrying two complete ORFs. Comparisons of DNA and amino acid sequences suggested the closely relationship between ERTBV and RTBV. The Oryza AA-genome species vary in the ERTBV copy number. The species carrying low-copy-number of ERTBV segments have been reported to be extremely susceptible to RTBV. The DNA methylation state of the ERTBV sequences was correlated with their copy number in the genome. CONCLUSIONS: These ERTBV segments are unlikely to have functional potential as a virus. However, these sequences facilitate to establish putative virus that provided information underlying virus integration and evolutionary relationship with existing virus. Comparison of ERTBV among the Oryza AA-genome species allowed us to speculate a possible role of endogenous virus segments against its related disease.

Characterization of the repetitive sequences in a 200-kb region around the rice waxy locus: diversity of transposable elements and presence of veiled repetitive sequences.

Repetitive genomic sequences might have various structural features and properties distinct from those of the known transposable elements (TE). Here, the content and properties of the repetitive sequences present in a 200-kb region around the rice waxy locus were analyzed using the available rice genomic database. In our previous Southern blotting analysis, 70% of the segments in this region showed smeared patterns, but according to the present database analysis, the proportion of repetitive sequences in this region was only 15%. The repetitive segments in this 200-kb region comprised 75 repetitive sequences that we classified into 46 subfamilies: 21 subfamilies were known TEs or repetitive sequences and 25 subfamilies consisted of newly identified TEs or novel types of repetitive sequences. The region contains no long terminal repeat (LTR) retrotransposable elements, but miniature inverted repeat transposable elements (MITEs) constituted a major class among the elements identified. These MITEs showed remarkable structural divergence: 12 elements were found to be new members of known MITE superfamilies, while five elements had novel terminal structures, and did not belong to any known TE families. Interestingly, about 10% of the repetitive sequences, including virus-like sequences did not have any of the usual characteristics of TEs, suggesting that a certain proportion of repetitive sequences that might not share the transpositional mechanisms of known elements are dispersed in the compact rice genome.